Hard sphere-like dynamics in a non hard sphere liquid 
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The collective dynamics of liquid Gallium close to the melting point has been studied using In- 
elastic X-ray Scattering to probe lengthscales smaller than the size of the first coordination shell. 
Although the structural properties of this partially covalent liquid strongly deviate from a simple 
hard-sphere model, the dynamics, as reflected in the quasi-elastic scattering, are beautifully de- 
scribed within the framework of the extended heat mode approximation of Enskog's kinetic theory, 
analytically derived for a hard spheres system. The present work demonstrates the applicability of 
Enskog's theory to non hard- sphere and non simple liquids. 

PACS numbers: 67.40.Fd, 05.20.Dd, 67.55.Jd, 61.10.Eq 



Since the first appearance of Enskog's theory [lj, spe- 
cial attention has been devoted to the theoretical and 
numerical study of the dynamics of hard spheres fluid 
as a useful tool to mimic the behavior of simple liquids 
During the seventies, the development of neutron 
scattering (INS) facilities provided a database of dynam- 
ical properties of simple liquids, while the enormous ad- 
vances of computational capabilities allowed recourse to 
the hard sphere model to evaluate transport coefficients 
and neutron scattering response ||. At the same time, 
on the theoretical side, the dynamical properties of an 
ensemble of hard spheres have been investigated by re- 
visiting Enskog's original kinetic theory. One of the most 
significant outcomes of these efforts is the so called ex- 
tended hydrodynamic theory [(J 0, H, @ , which has been 
particularly successful in describing INS data 01 ■ While 
Enskog's kinetic theory, strictly speaking, applies to hard 
spheres fluids only, it has the advantage of readily pre- 
dicting the Q-dependence of some transport coefficients, 
as opposed to more rigorous and involved theories using 
memory functions or generalized hydrodynamics. Thus 
this theory remains a very useful tool. 

The extent of validity of Enskog's kinetic approach 
has been satisfactorily tested against inelastic neutron 
scattering data collected in very simple liquids includ- 
ing Kr, Ar, Ne, Rb These are all similar in that 
their structure is well described by a hard sphere model, 
with an effective density and radius which can be de- 
termined by matching the first peak of the static struc- 
ture factor. However, no attempt has been made, to our 
knowledge, to verify the applicability of the theory to 
liquids markedly deviating from hard-sphere like struc- 
tures, probably because very few accurate, constant Q, 
experimental determination of coherent spectra exist in 
these systems [QJ. 

In this letter, we present a study of the collective dy- 
namic structure factor (S(Q,uj)) at wavevectors beyond 



the first maximum of the static structure factor (S(Q)), 
using inelastic x-ray scattering. The system we address 
is liquid gallium, as the purpose of this work is to ascer- 
tain validity of the connection between structural and dy- 
namical properties as predicted by Enskog's hard sphere 
theory for a liquid with strong non-hard- sphere struc- 
tural properties. Among simpler liquids, Ga has peculiar 
structural and electronic properties. In addition to the 
low melting temperature (T m = 303 K), it shows an ex- 
tended polymorphism in the solid phase with complex 
crystal structures characterized by the competition be- 
tween metallic and covalent bindings. Despite the nearly 
free electron electronic DOS, anomalies associated with 
some covalency residue have been reported [l2j |. Among 
them, the most important for the present purpose is that 
the first peak of the S(Q) presents a hump characteris- 
tic of non close-packed liquid structures 0, Q , which 
marks a significant departure from the smoother behav- 
ior of the hard sphere structure factor. Although the 
dynamics of li quid Gallium has been previously studied 
by INS mUllIl and IXS the restricted available 
Q range in the first case, and the presence of incoherent 
scattering in the second case, prevented a study of the 
collective dynamics in a Q region beyond the main peak 
of the structure factor like the present one. 

We demonstrate that the collective dynamics of molten 
gallium is well described by Enskog theory up to wavevec- 
tors as large as three times that of the first maximum of 
the static structure factor. While the reduced density 
and mean free path of liquid Gallium at melting point 
fall into the region of applicability of Enskog's theory, 
the structural peculiarities of this system would discour- 
age any recourse to predictions stemming from an hard 
sphere paradigm. However, the hard-sphere model turns 
out to describe the dynamics even in the region where the 
structure factor is notably different than a hard-sphere 
model. 
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The experiment was performed at the high resolution 
inelastic scattering beam-line (BL35XU) |l9j of SPring- 
8 in Hyogo prefecture, Japan. High resolution was ob- 
tained using the (9 9 9) reflection of perfect silicon crys- 
tals while a backscattering geometry (| — Ob ~ 0.3 mrad) 
was used (for both monochromator and analyzers) in or- 
der to obtain large angular acceptance. The flux onto the 
sample was « 10 10 photons/sec (100 mA electron beam 
current) in a 1.8 meV bandwidth at 17.793 keV. The use 
of 4 analyzers crystals, placed with 0.70 degree spacing 
on the 10 m two-theta arm (horizontal scattering plane), 
and 4 independent detectors, allowed collection of 4 mo- 
mentum transfers simultaneously. Slits in front of the 
analyzer crystals limited their acceptance to 0.24 nm" 1 
in the scattering plane. The over-all resolution of the 
spectrometer was about 2.8 meV, depending on the an- 
alyzer crystal. Typical data collection times were 200 
s/bin, where the bin size was fixed at 0.25 meV. 

The Ga sample, about 80 microns in thickness, was 
placed in sample cell with thin (2 x 250 micron) single 
crystal sapphire windows. This was held, in vacuum, at 
a constant temperature of 315 K. 

Measured IXS spectra are reported in Fig. ^ f° r se- 
lected constant values of momentum transfer Q. The 
strong central peak is readily visible along with some 
background due to phonons in the sapphire windows. 
Solid lines in the figure are best-fits with two components: 
i) a single Lorentzian line, modified in the usual way to 
account for detailed balance |20| , and convolved with the 
instrument resolution and it) the empty cell background 
measured at each Q, normalized by the sample transmis- 
sion. As it will be shown in the following, the Lorentzian 
spectrum is predicted by Enskog's kinetic theory, which 
also furnishes explicit expressions for the Q-dependence 
of the parameters describing its linewidth and amplitude 
(see, for instance, Eq. J^J). 

The main idea behind Enskog's theory is to evalu- 
ate the correlation functions of the microscopic quan- 
tities, such has the density-density correlation function 
of interest here, replacing the Liouville operator and 
the set of relevant dynamical variables defined at the 
N— particle ensemble level, with the one particle En- 
skog operator L and appropriate single particle vari- 
ables. Within this framework, the dynamic structure 
factor reads (|E0EI3) : 

with 

B j(Q) = (*? ) (Q,v 1 )) 1 <*f ) *(Q,v 1 ))i (2) 
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FIG. 1: IXS spectra of liquid Ga (T = 315 K) at the indicated 
fixed Q values (open dots). Also reported is the instrumental 
resolution (which is the same -dotted line- for all the reported 
spectra) and the best-fit lineshapes (continuous line, see text). 
The structures at larger energy transfers are phonons from the 
sapphire windows of the sample container. 

The subscript indicates a single particle averages 
over the Maxwell Boltzmann velocity distribution func- 
tion, while Zj(Q), ^ L) (Q,Vi) and ^.^(Q.vi) are the 
eigenvalues and the left and right eigenvectors of L, re- 
spectively. 

There are several approaches to determine the spec- 
trum of L, and different approximations can be applied 
depending on the density and kinematic regions of inter- 
est. Following Ref. , these regions are identified by the 
values of the reduced density 

V /V = P a 3 /V2 (3) 

in which Vq is the closest packing volume for spheres 
of radius a and number density p, and by Enskog's mean 
free path Ie = lo/x> with Iq = 1/ \/2pit(7 2 the Boltzmann 
mean free path and x = ff( cr ) the P arr distribution func- 
tion evaluated at the contact point between two spheres. 

Taking, as usual, density, momentum and energy as 
relevant variables, the lower three eigenvalues of L always 



3 



go to zero with Q — > 0. By introducing Enskog's thermal 
diffusion coefficient De, the adiabatic sound velocity c Q 
and the sound damping coefficient rjj,it can be shown 
that these low lying eigenvalues are 0, @ : 



zi(Q) = z h (Q) = D E Q 2 
Z2,a{Q) = z±(Q) = ±ic Q ■ 



r E Q 2 



(4) 



This limit is practically attained at low densities 
(Vq/V < 0.1 and therefore l E ~ lo) and sufficiently 
small Q's (Qa << 1). This condition normally occurs in 
the case of light scattering experiments from dilute gases 
(QIq « 1), and one speaks in terms of three extended 
hydro dynamic modes. 

In dense fluids {Vq/V approximately larger than 0.35), 
Kamgar-Parsi et al. have shown that a description in 
terms of three effective hydrodynamic modes still applies 
0, and the low Q limit of these modes is again coinci- 
dent with the hydrodynamic result. At variance with the 
previous case, however, the hydrodynamic scheme breaks 
down as soon as Ql E ~ 0.05. Above this value, Enskog's 
operator is dominated by binary collisions and only the 
extended heat mode is well separated by all the other 
modes 0. The following approximate expression for the 
extended heat mode can be given: 



in which De is Enskog diffusion coefficient and 



(5) 



d(Q) 



1 



l-j Q (Q(j) + 2j 2 {Q(j) 



(6) 



can be expressed in terms of the first two even order 
Bessel spherical functions. Enskog's diffusion coefficient 
is related to the Boltzmann diffusion coefficient 



of three Lorentzian lines (effective hydrodynamic modes) 
up to relatively large wavevectors (Q w 20 ran -1 ), 
while above this value the acoustic mode is overdamped 
and the quasi-elastic extended heat mode should domi- 
nate (QIe >> 0.05). The full width at half maximum 
(FWHM) of the quasi-elastic mode is quantitatively pre- 
dicted by Eq.s ©, JBJ and 0. 



> 



■4— > 

T3 




Momentum transfer [nm ] 



FIG. 2: FWHM of the dynamic structure factor (left axis): 
experimental determination (open circles, error bars smaller 
than the symbol size) and different predictions according to 
the revised Enskog's theory (Eq.JJ). Thin continuous line: 
a is the only fitting parameter (De is determined exploiting 
Eq.0). Thick continuous line: both a and De are fitting 
parameters. Thin dashed line (WA): a is determined as the 
value corresponding to the hard sphere structure factor which 
better describes the measured S(Q) (De is the only fitting pa- 
rameter). Thin dot-dashed line (MP): a is determined from 
the main peak of the structure factor, through Eq.® (De 
is the only fitting parameter). Also reported is the static 
structure factor (dotted line, right axis), which drives the os- 
cillations in the FWHM. 
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by the collision enhancing term g(a) as De — D Q /g(a). 
Plugging the well known analytic expression of g(cr) for 
hard spheres [2l| in the previous equation, one easily gets 
a final expression for De in terms of the packing fraction 

If = 7T/5(7 3 /6: 



D E = 



16 



TrksT 



6 (1 - V? 
m V (1 — y/2) 



(7) 



Summing up, by virtue of the results recalled sofar, 
the S(Q,u>) in dense fluids, such as in the present case 
of a monoatomic liquid close to the melting point, can 
be described to a sufficient extent of accuracy in terms 



The half width at half maximum of the quasi-elastic 
line (heat mode) are reported in Fig. 13 as a function 
of the Q values, together with the S(Q) [14j . Here, one 
can clearly observe oscillations in the line-width driven 
by the S(Q). The second minimum at Q ~ 30 nm -1 , in 
particular, is clearly related to the shoulder of the S(Q) 
observed the same Q position. The line-width predicted 
by Eq.© is also reported, with different choices of the 
two parameters a and De (see table QJ. By leaving both 
parameters free one gets a best fit value of a = 0.279 nm, 
corresponding to a packing fraction of tp = 0.601 (and a 
reduced density Vq/V = 0.811) which is just beyond the 
maximum theoretical value for hard spheres (ip — 0.545). 
This seems to be a reasonable result for a dense liquid, 
like the one under investigation here, recalling that a 
has to be regarded to as an effective parameter, thus not 
necessarily strictly related to the density through Eq. |J3J . 
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Using only a as free parameter (i.e. exploiting Eq.JTJ)) 
one still gets reasonable agreement (thin continuous line 
in FigEJ, even in the region of the secondary shoulder. 

Aiming at a description of microscopic properties of a 
simple fluid within the hard sphere paradigm, alternative 
choices of the effective particle diameter can be consid- 
ered. One can, indeed, choose a from the position of the 
first peak Qm of the structure factor (MP): 



or, through Eq.@, by adjusting the value of the re- 
duced density to maximize the whole S(Q) agreement 
[2^ | (WA). We tried both, either using Eq.0 or not, 
and the agreement with the experimental line-width is 
poor. The results are quantitatively summarized in table 
12 by fixing the hard- sphere diameter with one of the two 
mentioned criteria the \ 2 is definitely worse. Moreover, 
looking at Fig.J2J), one clearly observes that with this 
choices the secondary minimum observed in the FWHM 
at Q ss 30 nm _1 is smeared out. This feature, therefore, 
is not the mere consequence of the oscillations in S(Q) 
(De Gennes narrowing). It testifies to the presence of 
hard sphere dynamics and allows the sharp determina- 
tion of an effective hard sphere diameter via the Bessel 
terms of Eq.©. 

But what is the physical meaning of such "dynamical" 
effective diameter, and why it is larger than those ob- 
tained from the static structure factor? As we already 
mentioned, a description in terms of a uniform distribu- 
tion of hard spheres does not apply to liquid gallium, 
as well as to several other liquid metals (Zn, Cd, Bi, 
Si). The covalency residue of Ga, in particular, has been 
rationalized in terms of dimer-like structures 01 j an d 
clustering effects have been hypotized aiming to repro- 
duce the S(Q) The effective diameter that we find, 
therefore, could be an indication of the supra-atomic na- 
ture of the "effective particles". Consistently, plugging 
our effective a = 0.279 nm in Eq. Q, and solving it for 
a correspondent effective mass, we get m e ff = 83 a.m.u., 
which is larger than the atomic mass of gallium m = 69.7 
a.m.u. The effective mass and diameter, therefore, could 
be regarded as the mean sizes of the "clusters" undergo- 
ing collective dynamics. This hypothesis, however, needs 
to be tested extending the present study to other systems 
and at different temperature, as the structural anomalies 
are known to be strongly temperature dependent. 

In summary, we have shown how the microscopic dy- 
namics of a well known non hard-sphere liquid, namely 
liquid Gallium at the melting point, is surprisingly well 
described by Enskog's theory, analytically derived for a 
hard-sphere fluid. More specifically, the connection be- 
tween structure and dynamics, which is the main out- 
come of this theory, is robust enough to persist even in a 
wavelength region where anomalies typical of non simple 



liquids develop. From the experimental data, therefore, 
an effective hard-sphere diameter can be obtained, which 
turns out to be larger than the one related to structural 
properties according to well established methods. This 
result suggests, therefore, an extended validity of the re- 
vised Enskog's theory beyond the class of systems (hard 
spheres) for which it was derived, and provides an ex- 
perimental route to the determination of an effective dy- 
namical diameter. 
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TABLE I: Structural and dynamical properties of liquid gallium at the melting point, as derived by the present experiment. 
Each row is relative to a different fitting strategy, since the effective hard sphere diameter has been determined: i) as an 
adjustable parameter, ii) by the position of the main peak of the structure factor (MP), Hi) as the value corresponding to 
the hard sphere structure factor which better approximates the measured one (WA). Correspondingly, the Enskog's diffusion 
coefficient, De, can be obtained either as an additional adjustable parameters or fixed through Eq.Q. All the other quantities 
are derived by a and De according to the definitions given in the text. 



